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XYI. 



On Certain Forms of Interpolation, 



By W. P. G. BAETLETT, A. M. 



Communicated by Professor B. Peirce, May 13, 1862. 



The quantity to be interpolated being a function of t, the following scheme shows 
the notation and arrangement adopted for a series of its values and their differences. 
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Whenever the subscript numbers are omitted in writing differences, -|-^ is to be 
understood for all odd orders, and for all even ones. 

It may be observed that, if we put J_ i in place of Y , we shall have the following 
relations between its differences and the original ones of Y : 



(3) 



f 4 (<*_,) = 4, 4 (<*_,) =A, &c. ; 

and similarly, if we put J^ in place of F , 
4(4) =4, 4(4) = 4, &c. ; 



and so on with higher orders of differences, the law of this symbolic multiplication 
being to add the subscript numbers of the factors to form that of the product, as 
well as their (exponential) indices to form its index. 
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436 ON CERTAIN FORMS OF INTERPOLATION. 

Let i be a positive integer, and k = -; then after the insertion between Y n and 
Y n+X of » — 1 new values of the function (which is called interpolating to k ths ), 
let a portion of the series and its differences be 



The quantity & will be an implicit function of t, the 
values indicated in (4) corresponding to t = t — \ k, 

f o H~~ i ^' # o ~h f *> & c * » * et ^ ts va ^ ue f° r * = * i 1 be 
v N 5 ±i . 5^ will then, m aW cases, whether i is odd or even, be 

the same function of 5_j and its differences that 5 4 is of 

Y and its differences ; that is, if 



d±i = <P±t (^)> 

then 

»l = w (»-*) = w • »-* ( r o) = fo (To) ; 

( 5 ) -( and similarly for similar functions of higher orders : 

d" ±i = q> ±i (51) = <p ±i . qp t . <p_j (F ) = (j. 8 ± j (T ) , 

K = »i ( 5 -j) = <Pi • »-* • fi • f-i ( T o) = ft ( T o), &c. ; 

tf representing the result of the n successive functional operations indicated in each 
third member. 

From the well-known development 

™ _!_*„ ■ M*Tj) . , (* + !)* (A-l) . | (A + 1)A(A-1)(A=f2) 

(6) y h = y + A^ ± H 2 ^ "" 273 ±* "■ 2T3Ti ° + ' 

is directly obtained 

[l p 1 £2 32 ^2 "I 

" ±h ~T72Ts J> ±i + 17271 * 17175 ^±* ~ &C -J ' 

which is a linear function of certain differences, which are themselves linear functions ; 
therefore <p n is obtained by n symbolical multiplications of the series (7) according 
to the law (3). The subscript numbers are seen to result so that, omitting them from 
the operation, we have 5 n = the symbolical n th power of 8 ; that is 



(8) 



. r 1— v . . 1—e 9— p i 
= l l d u~ ' + -ai-- "so- '-J' 
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(9) , = v[s-l=**+l=*. *=**..], 

(10) * = *|/ _^» + -_. 24Q y...J, 

(11) * = *• [^ — g-^ + -J- . -Jg- J» ...J, 

(12) *• =p[j»-A(i_|?)^...]. 

It will be observed that the notation 8%? was substituted in (4) for each difference 
of the form.#^", because the symbolic even powers of 8 as there denned are always 
identical with €*&. When i is odd, 8 2n+1 is also identical with the middle iP n+l be- 
tween Y and Y 1 ; and all the series (8) ... (12) may be used directly. When i is 
even, let 8' 2n+1 denote half the sum of the two d 2n+1 s next above and below the 
middle point between Y and F^ ; then in each case it will be the same function of 
8 2n and its differences, that 8' is of Y and its differences ; that is, if 

(13) 3' =<p' (To), then 5' 2n+1 = <p< (6 2 n ). 
But we have for these functions of k, whether % i is odd or even, 

(14) 3' (h) = £ <$ (2 k) ; that is, <p' (k) = i <fi (2 k) ; 

therefore the series for 8' 2n+1 may be obtained by substituting 2 k for k in (7) and 
symbolically multiplying this into the series for J 8 2n . This gives 

7 T 1— -4P . . 1 — 4P 9 — 4F , 1 

(15) * =k [j 24— ' + -^4— * -80— ' ••]' 

r 1 — 2P . , /l — 2P 19 — 12 A* 1 \ . 1 

(17) 5l ^ ¥ y_t^ij...]. 

Besides the symmetrical relations obtained above by introducing the functions 8, 
the following may be obtained directly from (6) : — 

(18) ^ = 0j — eg, *_j = ©j + ©f; 

(# 4 + #_*= 0j + 0_j, 

in which 

, r i—* 8 , , i— a 2 2 a — p . „ i 

(20) e ±i = k [j ±i - -^ ^ ±i + —3- . -^-g- 4t j - &c.J , 

vol. viii. 57 
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, ni N » ,7/-. tn T o 1+* 2 — £ , , l+£ 2—£ 2-I-& 3 — £ . c 1 

(21) 0* = Ik (1 — i) [4 J_ . ___ 4 + _JT_ . ___ . _JT_ . ___ 4 _ &c j # 

Since the functions involved in (19) need not be restricted to integral values of i, 
similar expressions may be obtained for the differences between F and y ±m by sub- 
stituting m k fork in (9) and (20). By arranging these according to powers of mk 
and using their differences taken with respect to m, algebraic expressions may easily 
be got for any of the differences &l . To do this let m r be substituted for y m in the 
series (4) ; and, p? n m T being one of the n th order of differences of the new series, let 
&l (m) be the correspondingly placed difference of the series y m \ then 



(22) 



&l (— »») = *o" : ( m ) — ip o° ( m )> 

&% (m) = < ] ( m ) _j_ y/M ( w ) . 



in which, putting 

(23) jf-i = i (^ir 1 + 4"- 1 ) = s-r 1 + i 4 n , 

¥ ¥ 

*w (m) = r n m . k(j -i4 +A4S) + r m 3 .- (^ - i^) + r m 5 . — 4 + & c ., 

(24) j = f m . *4 - (^ - p- » 3 ) • y (4 - i4) ~ (^ ~ F"« 5 ) • ^o 4 + &c, 

*w (») = r m 2 • y (^ - tV4 + *V4) + r ™ 4 • 24 (4 - i 4) + r n ™° • 7^ 4 + &c 

Formula (22) is a perfectly general connection between the original and the interpolated 
differences of any order. 

When ra=l, (22) gives the means of deriving all the first differences,* # A , by 
computing the functions (24) for 5 or — — different values of m, according as i is 
even or odd, — in the latter case obtaining & = S^ from (8). For this method the 
second form of $ is adapted. But methods of this sort, however simplified in nu- 
merical application, seem unnecessarily laborious. 

Any interpolation may be made by combining (8) or (15) with (9), as follows. 
Compute the S 2 opposite each value of F, and the § or 8' opposite each interval. Then 
from the former fill in all the second differences, df, and by their means the first differ- 
ences, & k , and thence the function itself. If i is even, the two differences &-i and 

&+i, of which §' is half the sum, must be derived from 8' and the second difference, 

2 
#f, opposite it, by taking the sum and difference S'^^xH. This method is some- 

2 2 

times the most convenient, especially when i = 5. If the higher differences are 

* For another method of this sort, with a different system of differences, see Mk. Fekrel's article on 
Interpolation in the Mathematical Monthly, Vol. III. p. 377. 
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very large, the %P k s must themselves be interpolated, either by a distinct operation, 
or by the use of (10) or (16) and (11).* 

An objection to these processes, when many interpolated values are to be inserted, 
is that the sum of the computed first differences will not generally be exactly equal 
to the differences between two successive values of Y ; nor will the sum of the com- 
puted second differences exactly fill in the computed first differences, etc. These 
difficulties have to be remedied by subsequent arbitrary correction, or by com- 
puting the higher differences to one or more decimal places beyond those of Y, 
For this reason, interpolation to fourths, sixths, ninths, &c. should generally be 
done by successive interpolations to halves and thirds; and the choice of the 
special method of doing this must depend more or less upon the particular charac- 
ter of the function to be interpolated. The following are some practical methods 
with the numerical coefficients in the various series required. 

Interpolation to Halves. 

(21") el = l [4 - A A + A • A A • • . 

(15") d> =iJ, 

(9") « 2 =H^-tV^ + tV-£^---]- 

The interpolated differences may be obtained in successive pairs by taking the 
half-sum and half-difference of the quantities 

(s- t -\- &i = J = 28', 
(25) j it t 

When &{ can readily be obtained by interpolating the Sis, it may be more con- 
venient to use (9") instead of (21") to obtain it. 

Interpolation to Thirds. 

(«>'") «±* = *[^ ± | — A-^fci + A-'A^ti-"]. 

(21'") el = i [4 — 2 V 4 + A • t¥s 4 • • ; 

(8"<) 8 =\IJ —aV^ + A-i^---]» 

(9"0 « 2 =i[^--A-^ + A.A^...]. 

* This general method, as far as relates to odd values of i, with algebraic forms for 3j" -1 and dl n up to 
n = 4, is given by Encke in the Astronomische Nachriehten, Band 29, No. 695. Also a partial statement, 
obtained by induction apparently, of the symbolic connection of the series (8) . . . (12). 
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Formula (19) is sufficient to do the whole interpolation, using (20'") and (9'"). 
But it is often better (especially when, in interpolating to sixths or twelfths, it is 
desirable not to write down the inserted values of the function themselves, till the 
whole is finished, but only their differences) to use (8'") and (21'"), as follows : — 



(26) 



( &i + #i = A — d, &i = d, 

Interpolation to .Fourths. 

(20") ±i = i [J ±i - A J ±i + A . #A J ±i . . .] , 

(15") 3' =i|> -*-/ + *. A^---]. 

(9 IT ) « 2 =A[^-A^ + A.-ftV^...], 

(li") a* =rf»[^-A^-..]- 

Instead of interpolating twice to halves, & t and 5-, may be got by (19), using 
(20") and (9"), and &, and # f by 

U| + *j = 2«' = 4-(* 4 + * l ), 

(. ^J <^§ = ^2 • 

Either the second or third member of the first equation may be used; and &% 
may be got by interpolating 5 2 to halves, using for this purpose (11 IT ) if necessary. 

Interpolation to Fifths. 

(20 T ) e ±i = i [> ±i - A^t* + A • T 9 A4tj • • •] , 

(8 T ) 3 =*[^-A^ + A-t¥ 4 5^---L 

(9 T ) 5 2 = A l>" - A ^ + A • Air ^ • • •] , 

(io T ) a 3 = T a [^ - A ^ + A • tA ^ 7 • • •] , 

(ir) s* = jfal*—hJ>...]. 

If, besides getting ^ = 8 from (8 T ), we use (20 T ) and (9 V ) in (19) to get # t and 
#1, there will remain only #j and & t to be filled in. 



Q ±i 
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*^±i" 


_ 3 5 ,p 1 3 5 


143 


^ 5 ±i - 


■•]. 


$' 


= 


|[^- 2 V^ + 2V-i^- 


••]. 






«5 a 
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AC^- 


&** + *&■&&* 


...], 




<5' 3 
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-tV 7 ^ 5 +(tW 
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T^STs) 


z/ 7 . 
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Interpolation to Sixths. 

(20 TI ) 
(15-) 

(9 VI ) 

(16-) fl' s =, - 2 f e [yf - 1 ft ^ |- ( T ft . A - T -^ ro ,/ . . .] , 

(11-) 

If & t and ^ are computed by (19), and & t and & h by using (15 TI ) and (9 TI ), there 
will remain only # } and #| to be filled in. 

It will be observed that the series (8) . . . (12) and (15) . . . (17) always converge 
more rapidly than (20) and (21) ; and their use should sometimes be preferred on 
this account. A little reflection will often introduce, instead of the above-written 
fractions, some simpler equivalents, which may be mentally applied to the requisite 
degree of accuracy. Thus, in the examples given below, the following have oc- 
curred : — 

3 J. 1 7 1 1 



2T 



1 _1_ 1 14 I I 1 . 

— T T~ TOT j T2"5 — S H^ TT2T ' 



in all of which the right-hand term is easily applied when not small enough to be 
neglected. 

The function F, in the examples given below, is the Moon's declination for 1865, 
omitting the degrees ; the values are taken where their higher differences are quite 
large. Two methods of interpolating to sixths (by halves and thirds, and then by 
thirds and halves) are given, and their results may be compared. For computing a lu- 
nar ephemeris, the first method here given has been found to possess peculiar advan- 
tages. An example of interpolating to fifths, i. e. to every tenth of a day, is then given 
on Encke's plan. The whole of the work which it is necessary to write down is printed. 
The columns headed S . 8 &■, &c. contain the sums of each consecutive pair of values 
of the functions 8 Q 2 , &c. Since the sum and difference of any two numbers must 
be both odd or both even, in using (25) and (26) the nearest odd or even value has 
been taken for S . @ 2 , according as J and J — S were odd or even ; and similarly 
with 5 2 and &{ in using (19) and (9"). In the interpolation to fifths, the computed 
third and second differences, 6 s and d 2 , happen exactly to fill in the second and 
first respectively, all being carried to hundredths of a second, i. e. one place farther 
than the function Y. 



442 



ON CERTAIN FORMS OF INTERPOLATION. 



Date. 


Funct., Y. 


A 


A 2 


A 3 


J 1 


J 5 


Functions fob Interpolation to Halves by (25). 


A^ 4 


8 2 


8.8 e^ 


S.0 2 


d 

Apr. 27.5 
28.0 
28.5 
29.0 
29.5 


38 2.6 
58 56.8 
62 6.8 
48 20.1 
18 43.0 


1 II 

+20 54.2 
+ 3 10.0 
—13 46.7 
—29 37.1 


1 II 

—18 7.9 
17 44.2 
16 56.7 
15 50.4 

—14 31.4 


II 

+23.7 
47.5 
66.3 

+79.0 


+27.9 
23.8 
18.8 
12.7 

+ 7.4 


II 

—4.1 
5.0 
6.1 

—5.3 


+5.2 
4.5 
3.5 
2.4 

+1.4 


— 1093"l 
1068.7 
1020.2 

952.8 
— 872.8 


II 

—2161.8 
2088.9 
1973.0 

—1825.6 


II 

—270.2 
261.2 
246.7 

—228.1 





Functions 


foe Interpolation to Thirds by (19). 






Date. 


A S -^A* 


3 


-S.3 


tM 4 


9 5 2 


S.0 


S 2 


Apr. 27.5 
28.0 
28.5 
29.0 
29.5 


II 

+24.5 
48.5 
67.5 

+80.0 




i ii 
i-20 50.6 
- 3 2.8 
-13 56.7 
-29 49.0 


i n 

+23 53.4 
—10 53.9 
—43 45.7 


II 

+1.8 

1.4 

+0.9 


II 

—1066.0 

1018.1 

— 951.3 


i a 

+ 7 57.8 
— 3 38.0 
—14 35.2 


ft 

—118.4 

113.2 

—105.8 



Functions for Interpolation to Fifths by (8 v ), (9 t ), and (10 t ). 



Date. 



Apr. 



d 

27.5 
28.0 
28.5 
29.0 
29.5 



J-teS 



+24.2 
48.1 
67.0 

+79.6 



58 



--20 53.23 
--3 8.08 
—13 49.38 
—29 40.28 



+2.2 
1.9 
1.5 
1.0 

+0.6 



25 3 2 



-1090.1 

1066.1 

1018.2 

951.4 

- 872.0 



125 d s 



+24.2 
48.1 
67.0 

+79.6 



After Interpolation to Halves by (25). 


For Interpolation to every 2 Hours by (26). 


Date. 


A 


A 2 


A" 


A 1 


S 


S.9 2 


A-8 


S.0 2 


d h 

Apr. 27 12 

27 18 

28 
28 6 
28 12 

28 18 

29 
29 6 
29 12 


1 II 

+12 42.2 
8 12.0 

+ 3 45.6 

— 35.6 

4 50.0 

8 56.7 

12 54.5 

—16 42.6 


II 

—270.2 
266.4 
261.2 
254.4 
246.7 
237.8 

—228.1 


II 

+3.8 
5.2 
6.8 

7.7 

8.9 

+9.7 


n 

+1.6 
1.4 
1.6 
0.9 
1.2 
0.8 

+0.4 




1 II 

-2 44.0 
-1 15.1 
-0 12.0 

1 36.8 

2 59.0 
-4 18.3 


a 

—537.2 
528.2 
516.1 
501.5 
484.9 

—466.1 


1 II 

+5 28.0 

+2 30.5 

—0 23.6 

3 13.2 

5 57.7 
—8 36.2 


II 

—59.6 
58.7 
57.4 
55.8 
53.9 

—51.8 



ON CERTAIN FORMS OF INTERPOLATION. 



443 



After Interpolation to Thirds by (19). 


For Interpolation to every 2 
Hours by (9"). 


Date. 


& 


Function. 


J 


J 2 


z/ 8 


<5 2 


tf 


d h 

Apr. 27 20 

28 
4 
8 

12 
16 
20 

29 
29 4 


i a 
+4 58.1 
+2 59.7 

—0 52.4 
2 45.6 

6 24.7 
—8 10.5 


I n 

53 58.7 

58 56.8 

61 56.5 

62 59.2 
62 6.8 

59 21.2 

54 44.8 
48 20.1 
40 9.6 


i a 

+4 58.1 

2 59.7 

+1 2.7 

—0 52.4 

2 45.6 

4 36.4 

6 24.7 

—8 10.5 


II 

—118.4 
117.0 
115.1 
113.2 
110.8 
108.3 

—105.8 


II 

+1.4 
1.9 
1.9 
2.4 
2.5 

+2.5 


II 

—29.6 
29.3 
28.8 
28.3 
27.7 
27.1 

—26.4 


II 

—29.5 
29.1 
28.6 
28.0 
27.4 

—26.7 



Interpolation by using (25), and then (26). 


Interpolation by using (19), and then (9"). 


Date. 


& 


Function. 


& 


Function. 


d h 

Apr. 28 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 

29 


i a 
+1 44.6 

1 15.1 

45.9 
+0 16.9 
—0 12.0 

40.5 

1 8.7 

1 36.8 

2 4.5 
2 31.9 
2 59.0 

—3 25.8 ' 


58 56.8" 

60 41.4 

61 56.5 

62 42.4 
62 59.3 
62 47.3 
62 6.8 
60 58.1 

59 21.3 
57 16.8 
54 44.9 
51 45.9 
48 20.1 


1 II 

+1 44.6 

1 15.1 

45.9 

+0 16.8 

—0 11.9 

40.5 

1 8.8 

1 36.8 

2 4.5 
2 31.9 
2 59.0 

—3 25.7 


58 56.8 

60 41.4 

61 56.5 

62 42.4 
62 59.2 
62 47.3 
62 6.8 
60 58.0 

59 21.2 
57 16.7 
54 44.8 
51 45.8 
48 20.1 



Interpolation to Fifths by (10 t ), (9 T ), and (8 T ). 


Date. 


V? 3 


■a* 


& 


Function. 


d 

Apr. 27.9 
28.0 
28.1 
28.2 
28.3 
28.4 
28.5 
28.6 
28.7 
28.8 
28.9 
29.0 
29.1 


+0.27 
0.31 
0.35 
0.38* 
0.42 
0.45 
0.48 
0.51 
0.54* 
0.56 
0.58 

+0.60 


— 42".91 
42.64* 
42.33 
41.98 
41.60 
41.18 
40.73* 
40.25 
39.74 
39.20 
38.64 
38.06* 

—37.46 


1 II 

+2 44.57 

2 1.93 

1 19.60 

+0 37.62* 

—0 3.98 

45.16 

1 25.89 

2 6.14 

2 45.88* 

3 25.08 

4 3.72 
—4 41.78 


/ * ii 

56 12.2 
58 56.8 
60 58.7 
62 18.3 
62 55.9 
62 52.0 
62 6.8 
60 40.9 
58 34.8 
55 48.9 
52 23.8 
48 20.1 
43 38.3 



* These are the independent values derived from the original differences. 



